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We analyse the entanglement dynamics in a nonlinear Kerr-like coupler inter- 
acting with external environment. Whenever the reservoir is in a thermal vac- 
uum state the entanglement (measured by concurrence for a two-qubit system) 
exhibits regular oscillations of decreasing amplitude. In contrast, for thermal 
reservoirs we can observe dark periods in concurrence oscillations (which can 
be called a sudden death of the entanglement) and the entanglement rebuild 
(which can be named the sudden birth of entanglement). We show that these 
features can be observed when we deal with 2-qubit system as well as qubit- 
qutrit system. 
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1. Introduction 

One of the fundamental areas of interest in quantum information theory concerns problems 
of generation of entanglement in various quantum systems and the influence of external envi- 
ronment on disentanglement processes. As decoherence processes in entanglement evolution 
are unavoidable consequences of the system-environment interactions, one should consider 
this problem seriously in any real physical realisation of quantum computing devices. Es- 
pecially interesting in this context arc recently described phenomena of sudden death of 
entanglement [1, 2] which consist in entanglement disappearance in flnite time rather than 
exponential decoherence decay of an individual qubit, and a sudden birth of entanglement 
[3] characterised by the reappearance of entanglement after its death. Both these phenom- 
ena are caused by the interaction of the pair of entangled qubits with environment. For 
the systems described by Yu and Eberly [2] and also by Ficek and Tanas [3] an entan- 
gled pair was formed by the two 2-level atoms. The initially entangled atoms were located 
in separate cavities [2] (initially in the vacuum state) and had no direct interactions with 
each other during the evolution. A finite-time decoherence of the system having non-local 
properties is therefore caused by the spontaneous emission processes. Similarly, a sudden 
birth of entanglement in the system considered in [3] (composed of two initially entangled 
2-level atoms inside a single cavity in a vacuum state) is a consequence of interaction with 
the environment. The particularly, crucial role plays the formation of collective states of 
these two atoms followed by collective damping, which cannot be neglected whenever the 
atoms are close to each other. Moreover, a phenomenon of delayed birth of entanglement in 
the previously mentioned system was described in [4]. It was also shown that the thermal 
reservoir [5, 6] or a reservoir in a squeezed vacuum state [5] always causes sudden death of 
entanglement of a two-qubit system (two 2-level atoms) for certain classes of initial entangled 
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states. 

The decoherence process of entanglement between two harmonic oscillators coupled to 
the reservoir with the spectral density of the Ohmic character was considered in [7] . Various 
kinds of disentanglement evolution, depending on the characteristic features of that envi- 
ronment were considered. The both phases of evolution: (i) sudden death, (ii) sudden death 
and entanglement revival were identified as well as the usually observed phase of no sudden 
death behaviour. 

Some systems consisting of harmonic oscillators were also previously described in order 
to identify their various disentanglement properties. Within the Markovian [8] and for some 
cases under non-Markovian [9] approximations there is a possibility for obtaining a sudden 
death phenomenon, or for some cases in a non-Markovian reservoir the entanglement persists 
for longer times than in a Markovian reservoir. Under some conditions, the final state of 
the initially entangled pair can stay entangled (but with lower value of entanglement). It 
is also possible to obtain entanglement of the initially disentangled pair via the interaction 
with a Markovian reservoir [10]. 

In the present paper we will focus on the problem of evolution of entanglement generated 
between two anharmonic oscillators placed inside the two-mode cavity. We will show that 
it is possible to obtain a finite-time disentanglement (a sudden-death) of initially entangled 
qubits and for some cases also a revival of entanglement — a sudden-birth. We will consider 
not only the zero temperature reservoir but the thermal one as well, and we will show that 
the latter type of external environment is more suitable for observing the phenomena of 
sudden death of entanglement and its revival. 
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2. The model 

We deal with the model of Kerr coupler described in [11]. It was shown in that paper that 
when the initial state of the system was an excited one |2)a|0)b the system could be used 
as a generator of maximally entangled states (MES). From the quantum information theory 
point of view it was a qubit-qutrit system. The MES were generated most efficiently when 
the damping process of the generated photon states was negligible. The damping process 
usually causes decoherence in the system. This fact was shown on the basis of analysis of 
the time dependence of the fidelity between the quantum state generated and the Bell-like 
states. 

In the present considerations we will focus on the dynamics of the initially generated MES 
state (in the above described system) which is then left in a damped cavity. We will show that 
despite the damping processes the system is able to spontaneously generate entanglement 
after its sudden vanishing. In a two-atom system such a phenomenon is called a "sudden 
birth" of entanglement and it was discussed in [3]. A "sudden death" of the entanglement 
between two atoms was studied for example in [? ] . In our case the situation is qualitatively 
different as we deal with the entanglement between photon states and damping takes place 
within these photon states. 

The system under consideration is described by the following hamiltonian: 

Hnl + Hint + H^xt , (1) 

^{a^ra^ + fibTb^ , (2) 
e(at)262 + e*(5t)2^2 ^ 

aa^ + a*a . (4) 
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H = 



where 



Hnl — 

Hint = 
Hext — 



As usual, a(a^) and b{b^) are photon annihilation (creation) operators in modes a and b, 
respectively, Xa{Xb) are the nonlinearity parameters, a is the strength of the coupling between 
the external coherent field and the cavity mode a, and e describe the coupling between 
the two oscillators. The non-zero nonlinear coupling between the modes of the coupler is 
responsible for generation of the entanglement between the two-mode photon states. When 
the initial system's state is an excited two- mode one, for the values of e and a, small when 
compared with %, the form of hamiltonian (3) indicates that the states |2)(j|0)b and |0)a|2)b 
are involved in the system dynamics. Additionally, as a result of the action of the external 
field, the state |l)a|2);, must be also considered. As it has been proved in [11] the maximum 
entanglement between the states |2)a|0)6 and |0)a|2)6 can be formed and the entanglement 
between |2)(j|0)6 and |l)(j|2)ft states is also possible. 

3. The influence of damping 

As we want to describe the system in a damping environment we should use the density 
matrix approach. In the standard (Born and Markov) approximations, the adequate master 
equation can be expressed in the well known form: 




+ E [cljc,^ + c,jcl - clc,j - pc,^cl] , 

k=l 

where the operators Ck describe the damping in the modes a and b and are defined as: 



We have introduced the interaction with thermal baths of non-zero temperature via Ci^ and 
C'2„ terms, which include the mean photon numbers na and ni, of the system's environment, 
whereas Ci and C2 describe interaction with zero temperature bath. 

The degree of the entanglement can be described via various parameters. One of them is 
the concurrence defined for the two qubit system by Woofers [12] . The value of the concur- 
rence changes from for completely unentangled states, to 1 for the maximally entangled 
ones. This quantity is defined as: 

C{t) = max (yx[ - ^2 - - \/aI, 0^ . (7) 

The parameters Aj are the eigenvalues of the matrix R constructed from the density matrix 
p (obtained directly from the master equation) via the relation R = PcPc, where pc is the 
density matrix operator for 2-qubit system, and can be obtained from the relation: 

Pc^ Cry®CryK(^y®(Tv ■ (8) 

Appearing here Uy, is the well known 2x2 Pauli matrix. 

While for the two-qubit system the definition of the concurrence can be applied straight- 
forwardly, for a qubit-qutrit system there is a need to extract from the system's density 
matrix a matrix for a two-qubit system using a projection operator according to [13] . Thus, 

we get the reduced density operator p^ given by: 

Pc = no,2 o no,2pno,2 no,2 , (9) 

where 

no,2 ® no,2 = (|0) (0| + |2) (2|) ® (|0) (0| + |2) (2|) . (10) 
In such a way we will investigate the entanglement between the states for a two-qubit case 
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via the inspection of the formation of the two Bell-hke states: 

l^i) - ^(I2)j0), + i|0)j2)j , 

\B,) = ^(|2)jO),-z|0)j2)J . (11) 

One should remember that apart from these two states there is also a possibility of formation 
of other entangled states, which in general influence the total amount of entanglement in 
the system. 

We will analyse the process of sudden disentanglement and the possible reappearance 
of entanglement in two different situations. At first, we assume that the coupler initially 
prepared in an entangled state is located inside a cavity and we can ignore the external 
coherent field (a = 0). Then in the second case, we can allow the external field to interact 
with a coupler - in such a situation we supply the coupler with energy, which leaks through 
the cavity mirrors (one should remember that it is not the external field that produces the 
entanglement between the coupler parts). 

4. System without excitation (a = 0) 

Before we deal with the problem of entanglement in a qutrit-qubit system, first we will focus 
on the situation when the interaction with external coherent field is not present {a — 0). 
Physically, it corresponds to the faster (than in the o; 7^ case) leakage of the energy from 
a coupler system through the cavity mirrors, but moreover, it means that there are only 
a few states involved in the dynamics - for the case without damping there are the states 
|0)|2) and |2)|0), whereas, when the damping is present we also have to deal with these 
states:|0)|0), |0)|1), |0)|2), |1)|0), and |2)|0). 

When we include the damping processes, for the initially entangled state and for = 
nb = we can observe the continuous process of concurrence C {t) oscillations with decreasing 
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amplitude as expected - see Figl. We can identify this decoherence behaviour as a non- 
sudden death process and it is preserved for the vacuum reservoir case only. The whole 
population is finally transferred to the product state ;^|0)a <8) (|0) + 

For the thermal reservoirs we can observe that the decoherence changes in a different 
way (see Fig. 2) .We cannot observe the continuous and periodical vanish of the concurrence 
amplitude any more. The inset of Fig. 2 presents the part of the concurrence evolution for 
various values of Ua and Ub. Whenever in one of the bath modes (or in both of them) there is 
a different from zero mean number of thermal photons we can observe some gaps during the 
concurrence evolution. The concurrence may completely vanish — meaning that the states 
considered are not entangled any more — and after some time spontaneously appear again 
— the entanglement between the two states considered is rebuild due to the interaction 
between the Kerr-coupler and the thermal reservoir. 

Moreover, we have analysed the process of disappearance and revivals of entanglement for 
reservoirs with various values of Ua and Uf,. To simplify the numerical calculations, we have 
used greater value of the damping constants ■ja and 7b. For this case we can see that indeed, 
it is the interaction between the system and the thermal reservoir that causes a sudden 
disentanglement in the system, and moreover, it may cause a revival of the entanglement 
after that death. It can be seen from Fig. 3 showing maps of the concurrence for various 
numbers of photons. The plots are generated in such a way that dark areas correspond to 
C — value and the bright ones to any other (higher than zero) values of concurrence. The 
individual thin dark lines appearing in the maps reflect to the fact that C{t) during the 
oscillations reaches zero from time to time (in such a way as in Fig.l for example), and as 
such it is not interesting for our considerations. Only dark areas surviving for some time 
indicate that a sudden death has occurred. Any bright areas after that, indicate a revival 
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of the entanglement in the system. We can inspect this fact in more details analysing Fig.4 
showing the concurrence for Ub = 3 and for various Ua, with increasing values of Ua and Ub 
the time between the death and sudden birth of the entanglement becomes longer. 

5. System excited by an external field {a ^ 0) 

The second situation we want to focus on is that, when the external coherent field is switched 
on. This field's role is to preserve the energy in the system composed of two oscillators and 
the cavity, and to diminish the role of the leakage of energy through the cavity mirrors. One 
should remember that it is not the external field that causes the entanglement in the system. 
The entanglement arises due to the nonlinear interaction between the oscillators, which is 
present during the whole process [11]. An external coherent field apart from reducing the 
role of damping is also responsible for engaging other two-mode states in the dynamics. The 
entanglement in the coupler system without damping appears as a result of formation of the 
two Bell-hke states (11) and additionally, an entangled state of the form: 

I^3) = ^(|2)J0),+ |1)J2),) . (12) 

These states are formed alternately and the details can be found in [11]. When the coupler 
interacts with the reservoir, the whole population is finally transferred to the product state 
\P) = (|0) + |1))^ |0)b (in the time which depends on the damping introduced via the 
transparence of the cavity walls). Intermediate concurrence changes may appear, especially 
for the thermal reservoirs, and they may reveal previously mentioned features of sudden 
"death" and sudden "birth" of the entanglement, and in fact they can be even enhanced as 
the coupler is continuously supplied with the energy. 

As the reservoir is assumed to be zero temperature bath (in both modes a and b) we 
can observe the concurrence changes of a slightly different character. The two types of 
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concurrence maxima appear — see the inset oi at Fig. 5. These are two successive maxima of 
larger amphtude (they correspond to the formation of the BeU-hke states (11)) and between 
the two subsequent pairs of these maxima there are groups of maxima with significantly 
smaller amplitude (they are the effect of formation other entangled states — between |2)|0) 
and 1 1)1 2) states). In this way we deal with the physical system in which the previous 
entanglement is transferred out from the two-qubit system from time to time, because one 
of the states of the entangled pair |2) |0) is coupled to the remaining state 1 1) |2) of the qubit- 
qutrit system. One should keep in mind that the system interacts with the thermal bath 
that can significantly change the concurrence dynamics. This feature cannot be observed 
when the values of damping constants increase (see the inset (02) at Fig. 5). 

We have investigated the influence of the external fleld strength on the C{t) evolution for 
a thermal bath with 7^ and rib ^ 0. The exemplary map is presented in Fig. 6. We can 
see that even for a relatively small values of a there are explicit areas of no-entanglement 
and after some time the entanglement appears suddenly again. Moreover, in Fig.7 we can 
clearly see the effect of the thermal bath on the entanglement formation. For na = ni, = 1, 
when including the external coupling a, we can obtain a significantly enhanced period of 
disentanglement in the two-qubit system. The period of time when C{t) — become larger 
and larger with increasing values of a. For instance, when a = 0.3 (sohd fine in Fig.7) 
between t — 9 and t — lA the system is disentangled, and after t — lA the entanglement 
appears again. It means that during that period of time the entanglement is transferred to 
the state l^s). Next the population returns to the previously populated two-qubit system. 

6. Conclusions 

The possibility of changing the character of entanglement evolution caused by the external 
system's environment was considered. In particular, the entanglement obtained in a Kerr- 
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coupler system with nonlinear interaction between its parts, and the process of decoherence 
caused by a reservoir were analysed. At first the maximally entangled state in the system 
was generated, and after that the system started to interact with its external environment. 
The two cases, namely, that of the cavity interactions with environment of zero temperature 
{ua — rib — 0) a-nd the second one, when Ua and (or) rib are different from zero were 
discussed. The degree of entanglement was studied via the concurrence for an appropriately 
chosen two-qubit system. Moreover, analysis of the entanglement evolution for two physically 
different situations was performed. The first one was when the weak interaction with external 
coherent field was absent during the whole decoherence process and the second one, when 
such an excitation was present (for recollection — an external field does not generate an 
entanglement in a system). We have found for both cases that for a cavity in a vacuum 
state the process of entanglement decoherence is a continuous and manifested as oscillations 
of decreasing amplitude. Qualitative changes in the entanglement evolution for non-zero 
temperature external bath modes were observed. When there are additional photons inside 
the cavity, it was found possible to observe a sudden disappearance of entanglement in a 
two-qubit system, which can be called a "sudden death" of entanglement, and after some 
time its reappearance, which can be called a "sudden birth" of entanglement. 

When q; 7^ it is possible that, due to the interactions with photons inside a cavity, the 
entanglement can be transferred for some time to another entangled state \Bs) in which one 
of the two-qubit states is involved, and after some time, it returns to the system. 

For the case when a — there is no additional couphng to j^a) state. Different from 
zero number of photons in thermal bath leads to population of the other two-mode coupler 
states which (due to the interactions between the nonlinear oscillators) during the evolution 
change the population of the two-qubit system described. In such a way the interactions 
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with these additional two-mode states (and between them) lead to the population transfers 
from them to the Bell \Bi) and |i?2) states, causing the revival of the entanglement in the 
system after its disappearrance. 
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7. Figure Captions 

FIG. 1: The concurrence C{t) for the coupler initially prepared in a MES state and later 
left inside the cavity for Ua = rih = 0. Xa = Xb = X = 25, 'ja = lb = 0.001, a = and 
e = tt/x- The inset presents the same but in a different time scale. Time is scaled in 1/x 
units. 

FIG. 2: The concurrence C{t) for the coupler initially prepared in a MES state and later 
left inside the cavity for various values of na and and for 7^ = 7;, = 0.01. The remaining 
parameters are chosen the same as in Fig.l. The inset presents the same but in a different 
time scale. 

FIG. 3 The maps showing the dependence of C{t) on time and value of Ua for various values 
of rib. The other parameters are the same as in Fig.2. Black areas correspond to C{t) — 0, 
whereas the white regions to the positive values of C{t). 

FIG. 4: The concurrence C{t) for the coupler initially prepared in a MES state and later 
left inside the cavity for = 3 and various values of Ua- The remaining parameters are 
chosen the same as in Fig. 3. 

FIG. 5: The concurrence C{t) for the coupler initially prepared in a MES state and after 
that left inside the cavity for Ua — rib — 0. The parameters Xa — Xb — X — '^^j la — lb — 
0.001, a = tt/x ^ = ''^/X- The inset (oi) shows the same but in a different time scale, 
whereas (02) is plotted for 7^ = 7^ = 0.01. 

FIG. 6 The maps showing concurrence C{t) for various moments of time and values of 
external pumping a. The initial state is \Bi) , ria — rib — 1- The remaining parameters are 
chosen the same as in Fig.2. 

FIG. 7: The concurrence C{t) for the coupler initially prepared in a MES state. The other 
parameters are: Xa — Xb — X — 25, 7o = 76 = 0.01, e = tt/x and ria — rib — 1. 
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